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ABSTRACT 


The whole Sloan Digital Sky Survey (SDSS, 14, 555 deg^) has been searched 
for intergalactic globular clusters (IGCs) in the Local Group (LG). Using optical, 
infrared, and ultraviolet photometric selection criteria and photometric redshifts, 
the 2.1x10® of objects in the SDSS Galaxy Gatalogue were reduced to only 183,791 
brighter than tq = 19 that might be GGs. Visual examination of their SDSS 
images recovered 84% of the conhrmed GGs in M31 and M33 and yielded 17 new 
GG candidates, 5 of them of high conhdence, which we could conhrm as GGs in 
MegaPrime images from the Ganada, France, Hawaii Telescope. These 5 GGs 
are within M31’s halo, but the other 12 candidates are not close to LG galaxies 
or galaxies within 3 Mpc of the LG. Even though this search covers only one- 
third of the sky and some GGs could have been missed, it suggests that the LG 
does not contain a large population of IGGs more luminous than My ~ —6. In 
the direction of the M81 Group, the search yielded hve candidate GGs, probable 
members of that group. 

Subject headings: galaxies: groups (Local Group, M81) - galaxies: individual (M31) 

- globular clusters: general 
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Introduction 


This paper reports a search for globular clusters (GCs) over the whole footprint of the 
Sloan Digital Sky Survey (SDSS),14, 555 deg^ of the sky, out to the edge of the Local Group 
(LG) and beyond in the direction of the M81 Group of galaxies . It encompasses a large 


volume of the LG, beyond t 


Garrison-Kimmel et ah 


le virial radii of the Milky Way (MW) and M31 (~ 300 kpc, 


20141) . There are two reasons to suspect that GGs may be found far 


from the MW or M31 or any other LG group galaxy. These intergalactic globular clusters 
(IGGs) may have originated in galaxies, and through galaxy-galaxy interaction have become 
unbound from their ho sts and put on large orbits in the gravitational potential of the LG 


flBekki &: Yahagil 120061) . I GGs ma.y have also formed independently of galaxies in their own 


dark matter (DM) halos flPeebleslll984l) . While this second hypothesis has no compelling 
observational support at present, it has not been ruled out as a source of GGs in addition 


to the better supportec 

hvpc 

)theses that GGs formed c 

uring the early evolution of galaxies 

(e.g., 

Harris & Pndritz 

1994; 

Kravtsov & Gnedin 

20Q5; 

Rime green 

et ah 

2012 

) or from gas 

that was compressed w 

len galaxies collided (e.g., 

Schweizer 

1987; 

Ashman & Zepl 

1992 

J 


Whitmore fc Schweizer 


19951). 


IGGs hav e been discovere d in several ne arby galaxy cluste rs such a s Fornax flGregg et ah 


20091 ) . Virgo ( Lee et ah 


2010), Abell 118 5 ( West et ah 


2011 ). Goma ( Peng et ah 


2m, 


Abell 1689 (lAlamo-Martinez et al.l 120131 ). and in some of them the IGGs number in the 


thousands. The majority of the IGGs of these clusters have blue colors that are similar 
to the colors of many of the GGs found in dwarf galaxies and in the blue sequence of the 
typically bimodal color distribution of the GGs in massive early-type galaxies. The colors 
of these clusters are consistent with very old ages and low metal abundances. In the Goma 


cluster, however, a significant fraction (~ 20%) of the IGGs have been 
colors, suggesting that they are more metal-rich than the blue clusters 


bund to have red 


Peng et ah 


20111 ) . 
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The common interpretation of the presence of IGCs in galaxy clusters is that they formed 
in galaxies that were later either tidally disrupted or stripped of their outermost stars 
and GCs by galaxy-galaxy interaction. The galaxy merger Arp 105 in th e cluster Abell 


1185 , which is liberating both stars and GGs from their galaxy of origin flWest fc Gregg 


2ni4J ). is an interaction that is producing IGGs at the present time. The “hypervelocity 


cluster ”, which appears to have been ejected from the Virgo Gluster flGaldwell et ah 


may be further evidence of IGG producton by galaxy-galaxy interaction flSamsindl2015l) 


20 




Similar interactions ni ay have also occurred in galaxy groups according to the models of 
Bekki fc Yahagil (120061) . and the cluster GG-2 in the M 81 group, which l ies ~ 400 kpc from 


M81, may be an example of an IGG in a galaxy group flJang et ah 


2OI2I). 


Galaxy interactions have occurred in the LG, because the halos of both the MW 
and M31 show telltale signs of galaxy accretion. The ongoing disruption of the Sgr 
dwarf spheroidal (dSph) galaxy by the MW, and the substructures found in the halos of 
the MW and M31 provide the most direct evidence. Other evidence that the Galactic 
halo was formed by satellite accr etion has accumulated for more than 30 years (see 


Bland-Hawthorn fc Freeman 


of galaxy formation (e.g.. 


2014j. for a rev i ew). The simulations of the hierarchical pictur e 


Abadi et ah 


2003 


Bullock V .Tohn.ston 


2005 


Zolotov et ah 


200911 


show that the accretions of sub-halos of DM onto a large DM halo (i.e., a large galaxy) 
frequently lead to the capture of the sub-halo into a satellite orbit and eventually its tidal 
disruption. If the sub-halo is a dwarf galaxy containing stars and star clusters, then they 
become part of the stellar halo of the large galaxy. A few of the sub-halos that pass through 
the viral radius of the large halo are not captured but remain on orbits that are s everal time s 
the v iral radius in extent. Some of these so-called “back-splash galaxies” (e.g.. 


Gill et ah 


2005h may have had their gas removed by ram-pressure stripping and tidal stirring (e.g.. 


Mayer et af 


2001 


Kazantzidis et al 


20131 ) while passing close to the large halo. This may 


explain the existence of the Tucana and Getus dSph galaxies in the LG that are currently 
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far beyon d the virial rac 


their gas flTevssier et al. 


ii of M31 and the MW, and yet appear to haye been stripped o f 


20121) . The recent simulations by 


Garrison-Kimmel et ah! (120141) of 


galaxy groups containing two large halos that are similar to M31 and the MW indicate that 
more back-splash galaxies on large orbits are produced when two massive halos are present 
rather than just one. The case of the accretio n of a binary pair of sub-halos of unequal 


masses onto a large halo has been simulated by 


Sales et al. 


( 20071 ). who found that often the 


lower mass sub-halo was ejected at high velocity while the higher mass one merged with the 
large halo. It seems plausible then that during the accretion of a GC-bearing dwarf galaxy 
by a large galaxy, one or more of its GGs may become unbound from the dwarf and placed 


on large orbits, reminiscent of the orbits o 
for example the LG galaxies Fornax (e.g.. 


the back-splash galaxies. Some dwarf galaxies. 


Gole et al. 


20121 ) and NGG 6822 (Huxor et al. 
2013), contain GGs that are far from their centers and therefore may be easily stripped in a 


galaxy-galaxy interacti on. The large stellar substructures in the halo of M31 (iFardal et al 


2013 


Bate et ah! 


2014 and refs, therein) and it s large number of ha 


which appear associated with the substructures (jVeljanoski et al. 


o GGs, some of 


20141 ) suggest that it has 


of its halo GG system (e.g.. 


Zinn 


1993 


3 is true of the MW, as indicated by t 

le properties 

Mackey & Gilmore 

2004; 

Keller et al. 

2012 

) and 


rathe r directly by the accretion of GGs from the Sgr dSph galaxy (e.g.. 


Law fc Maiewski 


20101 ). It remains to be seen if the accretion events invo 
theoretical simulations predict were very numerous (e.g.. 


ving the MW and M31. which 


Bullock fc .Tohn.ston 


2005 ). placed 


any GGs on such large orbits that they are now IGGs in the LG. 

The hypo thesis that some GGs formed in their own DM halos is now three decades old 


([Peebles 


19841) . but it has not been confirmed or ruled out observationally. The discoveries 


of multiple stellar populations and variations in the abundances of He and other elements 
in GGs, has renewed the interest in this hypothesis as a way of explai ning how some GG s 


can retain gas and haye more than one episode of star formation (see 


Gonroy fc Spergel 
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201ll. and refs, ther ein). Because the DIV 


a large galaxy (e.g., 


Mashchenko fc Sills 


halo of a GC may be stripped as it is accreted by 


20051) . the most isolated and massiv e GCs should 


be the best candidates to prove 


1120111) and by 


his hypothesis. The recent investigations by 


Gonrov et ah 


Ibata et al.l 1 2013 ) of the structures of the GGs NGG 2419 and MGGl in 


the remote outer halos of the MW and M31 respectively, failed to hnd hrm evidence for 
the presence of DM. The 14 dwarf galaxies that are not members of the MW or the M31 


satellite systems (see 


McGonnachie 


20121) suggest the presence of isolated DM halos in the 


LG. In addition there might be DM haloes in the LG that formed single GGs. This idea 
and the possibility that some galaxy interactions may have flung GGs into larger orbits 
have motivated this search for I GGs in the LG, which great ly expands our earlier one of 


900 deg^ of the sky near M31 lldi Tullio Zinn fc Zinn 


2014 hereafter Paper II). 


To our knowledge, there has not been a previous survey of a large volume of the LG 
with the specihc goal of hnd ing IGGs. Many of the most rem ote GGs in the MW halo. 


including Pal 4 flAbell 


1955 ) and AM-1 ( Madore fc Arp 1979 ). both at distances from the 


Sun (c^q) > 100 kpc, were discovered by visual examination of the photographic plates 
taken with the Schmidt telescopes of the Palomar (POSS-I & -II), the European Southern 
(ESO), and the Siding Spring (U.K. Science Research Gouncil, SRG) Observatories. It is 
hard to judge to what distance these large surveys could have detected GGs because this 
depends on the properties of the clusters, the plate material, and the survey techniques. 
Our examination of the images of several of the brightest M31 GGs on him copies of the 
POSS-II plates and on images of the digital sky survey (DSS), suggests that these plates 
do not have the depth and resolution necessary to identify many GGs at th e distance of 
(Iq = 783 ± 25 kpc; all distances to LG galaxies in this paper are from 


M31 


McGonnachie 


1120121). A good exa mple is provided by the luminous cluster M31 MGGl (My = —9.2, 


Mackey et ah 


1120101 ). which was discovered on MegaGain ima. ges from the Ganada France 


Hawaii Telescope (GFHT) 3.6m telescope IjMartin et ah 


2006 ). On the him copies of the 
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POSS-II, MGCl appears stellar. On the SDSS images, which are the main source of our 
survey technique, there is no question that MGCl is a GC. Our survey can then identify 
similar clusters within the boundaries of the LG. Our survey is also sensitive to GCs, 
whose central regions are not resolved in the SDSS images, and are therefore catalogued as 
galaxies. On the other hand, it excludes GCs that are so resolved into stars that they are 
not classified as galaxies by the SDSS, and therefore do not enter our initi al selection of 


objec ts from its Galaxy Catalogue, such as the MW GCs Koposov 1 and 2 flKoposov et al 


20071) . Our survey misses also ultra-faint galaxies and dSph galaxies. Below we discuss the 


survey, its application, limitations, and results. 


2. Search Area 


Figure 1, which is a Hammer projection of the sky with equatorial coordinates, shows 
our search area (gray shading), the positions of the Galactic plane, the Galactic center, 
M31, M33, the other satellite galaxies of M31, the satellite galaxies of the MW, other 


LG galaxies, and galaxies with in 3 Mpc of t 


from the recent compilation of 


McConnachie 


l e LG. Most of these galaxies were selected 


(120121 ). which we augmented with the more 


2013b 

a) and the MW satellite, Crt I ( 

Belokurov et ah 

2 OI 4 I 


20141). which might be instead a 


low-luminosity GC in the outer halo of the MW flLaevens et ah 


20141) . 


In our first survey for IGCs in the LG, we surveyed ~ 900 deg^ of the SDSS in an area 
around M31 up to 500 kpc in projected distance from its center (i?gc) (see Fig. 1 in Paper 
II), which included most of its satellite galaxies. This area was selected because of the 
possibility that some of the accretion events that produced the large substructures and the 
rich population of GCs in M31’s halo, could have also produced IGCs. Our current survey 
(see Figure 1) encompasses larger expanses of the sky to the south and southwest of M31, 




























and it also expands in the opposite directions, with a large gap near the galactic plane. 
Note that the relatively low galactic latitude of M31 (—21?6) places it almost on the edge 
of the SDSS. Since the candidate backsplash galaxies Tucana and Cetus dSph galaxies he 
far from both M31 and the MW, IGCs may be similarly scattered throughout the LG. A 
wide distribution is also possible if IGGs formed in their own DM halos. For both reasons it 
is important to search as much of the LG as possible, and our current search covers about 
one-third of the sky. 


3. Survey Techniques 


The steps in this survev were modelec 

after our previous survevs for GGs in the 

remote halo of M31 ( 

di Tullio Zinn & ZinJ 

2013 

, Paper I) and for IGGs in the vicinity of 


M31 (Paper II). The hrst step in those two surveys was to select objects from the SDSS 
Galaxy Gatalogue, (a catalogue of non-point sources), with {g — i)o colors of old GGs: 

0.3 < {g — i)o < 1-5. The second step was to examine by eye the SDSS cutout images of 
these selected objects as provided by the SDSS website. The vast majority of the objects 
that passed the (g — i)o color cut could be immediately rejected as galaxies. Objects that 
were not easily rejected were then examined more closely in the r pass-band images, which 
we downloaded from the SDSS website in hts hie format. In the case of the IGG survey of 
Paper II, which covered a larger area than the survey in Paper I (~ 900 vs ~ 250 deg^), 
too many ambiguous objects still remained after this closer visual scrutiny. Gonsequently, 
additional steps were added in Paper II, which used combinations of optical, ultraviolet, 
and infrared colors to distinguish GGs from most galaxies on the basis of the shapes of 
their spectral energy distributions (SEDs). Although the survey in Paper II found several 
GGs, they are, however, more likely to be additional members of M31’s halo than IGGs 
because their Rgc values are < 140 kpc. The most time consuming part of these earlier 
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searches was the part that involved the initial visual inspection of the images (283,871 
images in the case of Paper II). Because the present search area is more than 15 times larger 
than the one in Paper II, it was impractical to employ the same initial technique of visual 
inspection. Instead, we were able to reject beforehand from the SDSS Galaxy Catalogue 
many thousands of galaxies by using some of the photometric selection techniques already 
experimented with in Paper II, based on their SED differences with GCs. The hnal step in 
the present survey was still a visual inspection of the SDSS cutout images, but of a greatly 
reduced number of objects. 


Since the halos of the MW and M31 are probably in part the debris from galaxy-galaxy 
interactions, the GCs in their halos may be representativ e of IGCs. The MW and M31 


GCs have similar colors and spectra fjSchiavon et ah 


2012 ). and the halo GCs in both M31 


and the MW appear to be older than 1 Gyr and metal-poor ([Fe/H]< —1). Consequently, 
for our search we chose color ranges that isolated similar objects. It is not clear that these 
broad ranges will necessarily include all IGCs that formed in their own DM halos because 
Peebles’ (1984) theory does not make specihc age and metallicity predictions beyond 
estimating that the formation process begins at a redshift ~ 50. We use the M31 GCs to 
serve as prototypes of IGCs, because their optical, infrared, and ultra violet colors can be 
measured with the same techniques that we will adopt for our survey. We will also examine 
how well our methodology identihes the conhrmed GCs in the other LG galaxies with GC 
systems. 

As we did in our two previous searches for GCs, the candidates were drawn from 


the SDSS Galaxy catalogue, which conta ins non-stellar o 


criterion. The SDSS frames pipeline (see 


jjects according to the following 


Stoughton et ah 


2 OO 2 I) measures the point-spread- 


function (psf) of each CCD frame and uses it to compute a “psf magnitude” for each 
photometric band, to which is applied an aperture correction to remove the effects of 
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variable seeing. The pipeline also measures the light-prohle of each object, which it hts with 
both exponential and de Vaucouleur prohles after they are convolved with the psf. The 
total magnitude of an object in each band is measured by the “composite model magnitude” 
(e.g., cmodelMag_r), which is based on a linear combination of the best htting exponential 
and de Vaucouleur prohles. The test for an extended source is made by summing over all 
bands the huxes captured by the psf magnitudes and separately the cmodel magnitudes, 
which are converted back to magnitudes. If the difference between these hux-summed 
magnitudes, psfMag—cmodelMag, is < 0.145, then the object is classihed as a staiS. 


Since all GCs will become indistinguishable from stars at some distance, it is important 
to investigate how the SDSS star-galaxy separation affects the completeness of our surveys. 
Figure 2 illustrates the effects of the star-galaxy separation and our magnitude limits on 
our surveys for IGCs in the LG and in the M81 group (see section 6). In each diagram, the 
open histogram is the luminosity distribution of the 168 GGs in M31, older than 1 Gyr, 
for which we could hud data in the “PhotoObjAH” catalog of the SDSS. These objects 


are conhrmed GG s according to versi on 5 of the Revised Bol ogna Gatalo, 


clusters (RBGv5, 


Galleti et ah 


( 20041) . 2012 edition) and/or 


? of M 31 globular 


Huxor et ah 


(12014 . and we 


have used the data from these sources and, for a few objects, from the SDSS to construct 
the histogram. In projected distance from M31 {Rgc), this sample spans 3.5 to 140 kpc, with 
80% at Rgc > 10 kpc. The solid histogram in the top diagram, which includes 146 GGs, 
87% of the total, are the ones that are classihed as galaxies by the SDSS and are brighter 
than our survey limit (ro = 19.0, see below). A sharp cutoff in My is not produced by this 
limit because some clusters of large radius and low central surface brightness (e.g, HEG8 
and HEGll) are measured systematically too faint by the SDSS. Because of this effect, 


there is a bias in our surveys against Ending similar ’’extended” clusters (see 


Huxor et ah 


^http: / / WWW. sdss.org/dr 12/algorithms/classify/#photo_class 
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2008|), which however appear to be rare objects (only 4% of the sample of 168 and 


of 


the conhrmed GCs in the RBCvS are classihed as extended). The ragged cutoff of the solid 
histogram is also caused by our tq limit depending only on the foreground MW extinctions 
of the objects, while the extinctions used to compute My include contributions from M31. 
In addition, 3 of the 168 clusters in the sample fail the star-galaxy separation and were 
listed as stars in the PhotoObjAll catalog. 


The solid histogram in the middle diagram, which contains 131 clusters or 78% of the 
total number, includes the effects of the tq limit and star-galaxy separation at 1100 kpc 
from the Sun. This may be approximately the limit of the LG in most directions of our 
survey because McGonnachie (2012) estimated that the zero velocity surface of the LG 
lies 1060 ± 70 kpc from the mid point between the MW and M31. To model the SDSS 
star-galaxy separation, we used the composite profile £t in each band to find the radius of 
the aperture that enclosed the same flux as the psfMag. This radius was typically 80% 
of the full-width-half-maximum (FWHM) of the psf as listed in PhotoObjAll. We then 
scaled the effective radii of the exponential and de Vaucouleurs profiles that are listed in 
the PhotoObjAll catalog from the M31 distance to 1100 kpc and computed new composite 
profiles for each band. Using these profiles and the effective apertures of the psfs, we 
computed the difference psfMag—cmodelMag, which, if < 0.145, indicated a stellar source. 
According to this procedure, 10 of the 168 clusters would be considered as stars at a 
distance of 1100 kpc, but only 2 of them are brighter than the rg cutoff. 


The solid histogram in the bottom diagram of Figure 2 shows the effects of the rg 
cutoff, now increased to 20.0 (see section 6), and the star-galaxy separation at the distance 
of the M81 Group, 3.6 Mpc (Karachentsev et ah 2013). Only 35 clusters or 21% of the 
original 168 are included in this histogram. Nonetheless, it suggests that ~ 50% of the most 
luminous GGs (My < —7.8) are non-stellar at M81’s distance. This result is in qualitative 
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agreement 
studied by 


with the presence of some M81 GCs in the SDSS galaxy catalogue (e.g., the two 


Jang et ah 


(I2OI2I)). 


For the optical data, we downloaded from the Galaxy catalog the SDSS ’’model 
magnitudes”, which are based on the better one of the exponential or the de Vaucouleur 
fits to the r-band light profile. While the model magnitudes provide the best measurements 
of the colors of galaxies, they are less good than the cmodel magnitudes for measuring 
the total lighio. Over the magnitude ranges of interest here the differences between these 
magnitudes are so small that they can be safely ignored (e.g., for the 147 M31 clusters with 
To < 19.0 in Fig. 2, the difference, modelMag—cmodelMag, have means of 0.031 ± 0.010, 
0.016 ± 0.005, and 0.008 ± 0.009 in the g, r, and i bands). We only selected objects with 
reddening corrected r magnitudes 12.0 < tq < 19.0, because our experience with the M31 
GGs suggested that fainter ones might fail our color selection criteria, and with very few 
exceptions, cannot be conhrmed as star clusters by visual inspection of the SDSS images. 
For an optical color we chose g — i because it is well measured for these relatively bright 
objects and has a long wavelength baseline. In addition this time we selected for our sample 
objects with reddening corrected colors in a more restricted range, 0.3 < {g — i)o < 1.1, 
which still encompasses the range in color of the outer halo GGs in M31 and most of the 
GGs belonging to the dwarf galaxies of the LG (see section 4). Because of the age-metallicity 
degeneracy effect, this range in color encompasses wide enough ranges in both metallicit y 


and age,—2 .25 < [Fe/H] < —0.33 if age > 1 Gyr, according to the models of 


(11 9981 . 


Marasto: 


]■ 


20051 ). Moreover the selection criteria for GGs that we developed in Paper II, using 


color combinations from the SDSS, the Wide-held Infrared Survey Explorer (WISE), and 
the Galaxy Evolution Explorer (GALEX) satellites, work well within this range of {g — i)o, 
but not outside it, as shown in Figure 3. These criteria were presented and discussed at 


2 WWW. sdss. org/ dr 12 / algorithms / magnit udes / 
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greater length in Paper II. Because the recent Data Release 10 (DRIO) of the SDSS Galaxy 
Catalogue is linked to the WISE All-Sky Catalogue, much of our present survey could now 
be done automatically via Casjob on the SDSS website. 


Our previous experience in searching for GCs in the halo of M31 has shown that the 
largest contamination in our sample comes from galaxies resembling GCs in appearance in 
the SDSS images because they look compact and nearly round. When the spectroscopically 
measured redshifts (z) of these galaxies were available, they indicated that the galaxies lie 
far beyond the LG. In most cases, the SEDs of these galaxies do not closely resemble those 
of GCs. Unlike GCs, many of them are still forming stars. Ones that have little or no star 
formation have stellar populations that are more metal-rich, and hence redder than the 
GCs found in the halos of the MW and M31 or in dwarf galaxies. Therefore these groups of 
galaxies can be separated from GCs by means of photometric criteria. To illustrate how our 
selection criteria can distinguish between these galaxies and GCs, we have selected a sample 
of 122 GCs in M31 with Rg^ > 10 kpc and 12 < tq < 19, that do not lie in close proximity 


to bright stars or in the very dense star he 


selected from the catalogues of 


Kang et ah 


ds o: 


20121 ); 


M31. These c onhrmed cl usters were 
( 2008 . 20 iJ l and Papers 


Huxor et ah 


I and II. Because we are searching for IGCs, in the construction of the selection criteria 
we selected a more remote sample {Rgc > 10 kpc) than the one used in Paper II {Rgc > 3 
kpc). The clusters in the inner regions are more likely to have formed within M31 rather 
than in dwarf galaxies or in their own dark matter halos. In the diagrams of Figures 3 and 
4, this sample of GCs is compared to the objects in a small area of our survey covering 
4deg^, and centered at RA, Dec = 161.0, -1-1.0 (1 = -1-249.1, b = -1-50.1). This region is 
representative of all but the lower galactic latitude regions of our survey, where the SDSS 
Galaxy Catalogue contains fewer galaxies per square degree and many more tight groups of 
stars. In this comparison region, the SDSS Galaxy Catalogue contains 2146 objects with 
12.0 < tq < 19.0. The SDSS cutout image of each of these objects was visually examined 
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and none of them was a candidate GC according to the criteria that we discnssed in Paper 
I. We nsed these objects to illustrate how the majority of similar objects can be separated 
from GCs without resorting to a preliminary visual inspection, which would have been very 
time consuming for the area of 14, 555 deg^ surveyed in this search. 


The top diagram of Figure 3 shows that many galaxies (X’s) are redder in the {g — i)o 
and {i — Wl)o colors than GGs (circles). The color {i — Wl)o is formed from the magnitudes 
in the SDSS i-band, transformed to the Vega system, and the WISE W1 band, which is 
centered at 3.4/im, while the W2 and W3 bands are centered at 4.6 and 12/xm. We used 
the prohle htted magnitudes for all WISE measurements because they are presented as the 
most reliable mea surement for unresolved sources in the description of the WISE catalogue 


(iGutri et ah 


2011 ). Nearly all of the M31 GGs are unresolved by WISE in the Wl, W2, and 


W3 wavelength bands, which have psfs with full-width-half-maxima (FWHM) of 6"1, 6"8, 
and 7"4, respectively. Later we will discuss the effects on our survey of forming colors with 
the SDSS model magnitudes, which include nearly all of the light of the objects, and the 
WISE and GALEX measurements, which do not. The top diagram in Figure 3 shows, 
in the {g — i)o — (i — IFl)o plane, the locations of the objects of the comparison region 
(X’s), which extend beyond the limits of this plot, and of the selected sample of M31 
GGs (circles). The dashed contour in this diagram encloses 90% of the density of the 2-D 
kernel density estimate of the GGs. The rectangular area, which encloses 120 (98%) of the 
GGs, marks some of the cuts that we will impose to isolate GGs, 0.3 < {g — i)o < IT and 
(i — IFl)o < 2.03. Only 389 (18%) of the objects in the comparison region pass these cuts. 


The middle diagram in Figure 3 shows the usefulness of adding a second cut with the 
photometric redshift, Photoz, which is derived from the 5 SDSS photometric bands by 
the K-D tree method and listed in DRIO. In Paper II, we showed that Photoz is a useful 
discriminant between M31 GGs and galaxies because it systematically assigns low values 





15 


of z to the GCs (see Fig. 4 in Paper II). Although these Photoz values are much larger 
than the real z’s of the GCs, they are smaller than the ones of many galaxies of similar 
color and magnitude. For this study, we chose for Photoz a wider limit (< 0.13) than in 
Paper II to reduce a bias against fainter GCs, like B407 which would not pass otherwise 
our selection. The original sample of 122 M31 GCs had to be reduced to 89 in the middle 
diagram of Figure 3, because values of Photoz were available for only the ones within the 
SDSS footprint. One of these clusters, B517, has an exceptionally large value of Photoz, 
0.177, and it was excluded when calculating the 2-D kernel density and its 90% contour. 
The {u — g)Q color of B517 is approximately 0.2 mag smaller than the ones of other M31 
GCs of similar [g — i)o, which may be responsible for its anomalously large Photoz. The 
rectangle in this diagram indicates the region delineated by the cuts in Photoz and {g — i)o 
that we will use to isolate GCs, and it encloses every GC (circles) except B517. The 389 
objects in the comparison region that passed the cuts in the upper diagram are the ones 
plotted in this middle diagram; 338 of them passed the Photoz cut. 

The bottom diagram in Figure 3 is a plot of the color W2-W3, versus {g — i)o. W2-W3 
is sensitive to the amount of star formation because the W3 band includes the emission 
from polycyclic aromatic hydrocarbon molecules. While all 122 M31 GCs are plotted in 
this diagram, W3 was measured for only 31 (solid circles). For the other 91 GCs, W3 is 
an upper limit (open circles), and therefore W2-W3 is also an upper limit. Using only the 
GCs with measured values, we calculated the 2-D kernel density, and the dashed contour 
in the diagram encloses 90% of the density of these GCs. The faintest GCs are the ones 
lying above the contour. In order not to exclude faint GCs with only upper limit values, 
we will impose two different cuts by W2-W3. If W2-W3 was a measured value, like for 
solid circle objects, the cut will be at 3.1. If W2-W3 was an upper limit, it is raised to 4.0. 
The majority of M31 GCs (~ 96%) passed all of these three criteria. Also plotted (crosses: 
upper limit; X’s: measured value of W2-W3) are the 338 comparison objects that passed 
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the previous two criteria. Only 128 of them or 6% of the original 2146 objects passed 
the hnal cut, demonstrating the value of this sequence of selection photometric criteria to 
reduce galaxy contamination. 


This selection procedure through WISE can be largely automated because the SDSS 
Galaxy Catalogue is linked to the WISE All-Sky Catalogue. Therefore, the first step for 
our selected survey for GC in the whole SDSS footprint was to query the SDSS Galaxy 
and the WISE All-Sky catalogues, via CasJobs on the SDSS website, for objects with 
12.0 < ro < 19.0, 0.3 < (^ - i)o < 1.1, {i - Wl)o < 2.03, Photoz < 0.13, andiy2 - 1^3 < 3.1 
for measurements of W2-W3, or < 4.0 for upper limits. Since the footprint of the SDSS 
contains ~ 2.1x10^ galaxies of all magnitudes and colors, a huge number of objects 
4.21x10^) passed our first five selection criteria, still a sample too large to be inspected 
visually in a reasonable time, and still containing a large component of galaxies. At this 
point we introduced another cut, also discussed in Paper II, based on the color [NUV — (/jo, 
which demonstrated to be very useful for separating GCs from galaxies with large UV 
fluxes, presumably because they contain many hot and luminous young stars. NUV is 
the GALEX magnitude in its near UV band (A effective = 2271 A) and g is the SDSS 
g-band magnitude. This additional criterion is illustrated in Figure 4 with the same 
sample of M31 GCs previously used (open circles) and the objects in the comparison region 
(crosses and X’s) that passed the previous cuts, and could also be identihed in the Galex 
ca talogue. Th e NUV measurements for many of the GCs were taken from the catalogue 


of 


Kang et al. 


(1201211 . For the remaining clusters and for the whole samples of comparison 


galaxies, we compared their positions with the objects in the GALEX catalogue that is 
available through the GALEX view website and considered a match if their positions agreed 
to < 6'.'0. The average difference between the SDSS and GALEX positions was only T.'4, 


with < 10% having differences > 3". Conseque ntly, they are srn a 


cases, comparable to the NUV psf (5'.'3 to ~ 8", 


Morrissey et al. 


ler t 


ran or, in the worse 


(120071) ). These results and 
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the relatively low density of GALEX sources on the sky (e.g., 


Bianchi et ah 


20141) . suggest 


that very few of the objects drawn from the SDSS catalog were matched with the wrong 
GALEX source. If more than one GALEX measurement was available for an object, we 
computed the weighted average, using as weights one over the square of the error listed for 
NUV. Using the sample of 122 GGs, we computed the 2-D kernel density and the dashed 
contour in Figure 4 encompasses 90% of the density. The curve that separates most of the 
comparison objects from the GGs is identical to the one plotted in Figure 2 of Paper II, 
which had a larger sample of galaxies. While all 122 GGs lie below this curve and therefore 
pass this last criterion, only 16 of the remaining comparison objects lie below, and none 
of them could be classified as candidate GGs by visual inspection. After these tests, we 
concluded that the combination of the selection criteria illustrated in Figures 3 and 4 could 
reduce the number of objects to be visually examined to a small percentage of the initial 
sample (~ 1%), while retaining most of the GGs (96%), making feasible our planned survey 
of the whole SDSS footprint. Therefore the second step, based on the color {NUV — g)o of 
our procedure was to query the GALEX view website for GALEX sources that matched the 
positions of the objects (~ 4.21x10®) that passed the WISE and Photoz cuts introduced in 
the first step of the survey. We found approximately 3.17x10® coincidences, and to these 
objects we applied the cut based on the {NUV — g)o color (Fig. 4). The objects that passed 
this last cut, added to the ones that we could not hnd in GALEX, formed our hnal sample 
of 183,791 objects to be visually inspected. 


4. Application to the Globular Clusters in other Local Group Galaxies 

The above criteria appear to be well tuned for separating the halo GGs of M31 from 
galaxies, but do they also work for the GGs in the smaller galaxies of the LG, which may 
be better analogues of the IGGs? To examine this question, we compiled similar data 
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on the GCs in the Fornax dSph galaxy, a satellite of the MW, in fonr satellites of M31 
(M33, NGC 205, 147, 185), and in two dwarf irregular galaxies of the LG, NGG 6822 and 
Wolf-Lundmark-Melotte (WLM). Because most of these galaxies do not lie in the SDSS 
footprint, we could not include measurements of Photoz in this discussion. Moreover we did 
not consider the GGs in the MW, the Magellanic Glouds, and the Sagittarius dSph galaxy 
because they subtend such large angles on the sky that their measurement would require 
substantially different techniques than we employ in our survey of the LG. 


In the case of the Fornax dSph galaxy (c^q = 147 ± 12 kpc), its GGs subtend moderately 
large angles on the sky, w ith Rh values of 477, 674, 770,1079, a nd 1877 for GGs F4, F3, F5, 


F2, and FI, respectively (jMcLaughlin fc van der Marel 


20051) . They provide a good test 


of the effects of using magnitudes, such as the WISE and GALEX ones, which do not 
capture all the light of an object. We computed the {i — lFl)o and W2-W3 colors using 
both the profile-htting WISE magnitudes, which we adopted for our survey, and the largest 
WISE aperture, which is 24775 in diameter. For the Fornax GGs, we could not hnd in 
the literature measurements in the g or i pass-bands. Gonsequently, we made estimates 


that a re based on the V and B-V measur ements of F2, F3 , F4, a nd F5 by 


Gordon fc Kron 


1119831). on the valu e s com piled for FI by 


Harris fc Racinel (119791) and the reddening from 


Mackey fc Gilmorel 1120031 ). We then constructed a transformation between {B — V)o and 


[g — z)o from measurements of M31 GGs in both colors. 


In the diagrams of Figure 5 are plotted the Fornax GGs and the GGs in the other LG 
galaxies of our sample. The contours are the ones calculated for Figures 3 and 4. For each 
Fornax cluster, the colors formed by the profile-fitted WISE mag (hlled circles) and the 
largest WISE aperture (X’s) are plotted and joined by dashed lines. Given the large Rh of 
FI, it is not surprising that with the profile-htted W1 magnitude, {i — W1 )q lies far outside 
the 90% contours that are dehned by the sample of M31 GGs. FI nonetheless passes our 















19 


{i — Wl)o criterion. Despite the fact that the Rh of FI is signihcantly larger than the psfs of 
W2 and W3, its W2-W3 color passes the cnt when the prohle htting magnitudes are used. 
The results obtained for other Fornax GCs with prohle-htting magnitudes are less deviant 
than the ones for FI, and with the largest WISE apertures all 5 Fornax clusters lie within 
or very close to 90% contours of the M31 sample. In the plot of {NUV — g)o vs. {g — i)o, 
all 5 Fornax clusters lie within the 90% contour despite their relatively large angular sizes. 
Because Fornax is well within the boundary of the MW’s halo as defined by its satellite 
galaxies (see McConnachie 2012), FI may be larger in angular size than any IGC in the 
LG, unless they are atypical GGs. None of the GGs in the more distant LG galaxies have 
Rh > 10" (see refs, below). The test provided by the Fornax clusters in Figure 5 suggests 
that the mismatches in the amounts of light included in the SDSS model magnitudes and 
the WISE and GALEX magnitudes are unlikely to cause GGs to be missed by our selection 
criteria. 


Recent surveys of t he dwarf irregular g alaxy NGG 6822 (cIq = 459 ± 17 kpc) by 


Hwang et ah 


(2011) and 


Hnxor et ah 


(120131 ) have added 7 new GGs, bringing the total 


known to 9. We could reliably place 6 of these clusters (HVII, G2, G3, G4, SG6, SG7) in 
the diagrams that use WISE measurements and 4 of them (G2, G3, G4, SG7) in the one 
that uses GALEX (solid triangles in Fig. 5). The other clusters are too contaminated by a 
rich star held (HVIII), a nearby bright star (Gl), or simply too faint (SG5) to be measured 
by WISE and GALEX. The identihcation of HVII and SG6 in the GALEX catalog was 
problematic because of other nearby sources. The above references provided the g and 
i-band photometry and the reddening values that we used. For HVII, we estimated (g — i)o 
from the (V — I)q color that was listed by Hwang et al. (2011). The reddening is variable in 


the direction of NGG 6822. 


Hwang et al 


(120111 ) quote yalues of E{q — i) 


for the clusters in their study. The {g — i)^ color of G3 (IHwang et 


al. 


red (1.35) for its metal-poor composition([Fe/H] = -1.61), which 


rom 0.16 to 0.36 


201111 is exceptio nally 


Hwang et al 


(l2014h 
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measured from its spectrum. C3 failed every one of our GC criteria because of its red color. 
Our survey technique would have discovered 5 of the 9 clusters in NGC 6822, (or 56%). 
This limited success should not be indicative of our search for IGGs in the LG because some 
of the failures and omissions of the NGG 6822 clusters were caused by contamination from 
other sources in this galaxy, and presumably IGGs are in less dense helds. 


The other LG dwarf irregular galaxy with GGs is WLM {(Iq = 933 ± 34 kpc), which 


contains one 
measured by 


uminous GC {My ^ —8.9). We estimated its {g — i)o from the (V — I)o color 


Billett et ah 


(120021) . It is identihed unambiguously in the WISE catalog, and 


passed our criteria based on WISE measurements (open circle in Fig. 5). It was not found 
by GALEX, presumably because of its proximity to other UV sources. 


M33 ((Iq = 809 ± 22 kpc), the spiral galaxy companion of M31, lies within the SDSS 


footprint, although few of its GGs are in the SDSS Galaxy catalog because they are proiec 


ted 


on the dense star helds of t 


(I2007h : 


Huxor et al. 


(I2009fl : 


lis face-on galaxy. Usi ng the lists of 


Cockcroft et al. 


Saraiedini &: Mancone 


(201 ll ). we identihed 10 GGs in the SDSS 


Galaxy catalog with reliable photometry, 7 of which are brighter than the vq cutoh of our 
survey and passed als o our Photoz criterion . All 7(crosses in Fig. 5) passed the (i — Wl)o, 


and all but one (S in 


Cockcroft et ahl 


201 ih . passed the W2-W3 one. Five of the seven 


clusters were identihed in the GALEX catalogue, and all of them passed the {NUV — g)o 
criterion. 


Veljanoski et al. 


(2013) recently discovered new GGs in the M31 satellites NGC 147 


(d© = 676 ± 28 kpc) and NGC 185 (d© = 617 ± 26 kpc) and provided g- and i-band 
photometry, which we used for both the new and the previously known GGs. Of the 10 
GGs now known in NGC 147, one (SD-GCIO) is fainter than our tq limit and another 
(PA-N147-3) had to be excluded from our analysis because the g and i measurements for 
it are only lower limits (see Veljanosky et al. 2013). Five were identihed unambiguously in 
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the WISE catalog (Hodge I & III, SD-GC7, PA-N147-1 & 2) Two of the remaining GCs 
(Hodge H & SD-GG5) appear to be confused with other sources, and the faint cluster 
Hodge IV was not recognized by WISE, perhaps because of the high stellar density in its 
field (see Veljanosky et ah 2013) . All 5 measured clusters (solid squares in Fig. 5) passed 
the {i — Wl)o and W2-W3 criteria. Four of them were identihed in the GALEX catalog, 
and they all passed the {NUV — g)o criterion. 


In NGG 185 there are now 8 known GGs f|Velianoski et al 


20131) . and 6 (FJJ-I, HI, IV, 


V, VH, PA-N185) were found in the WISE catalogue. FJJ-H lies in a dense stellar held, 
and FJJ-V HI appears to be blended with a nearby source. All of the 6 clusters with WISE 
measurements (open triangles in Fig. 5) passed the (i — iyi)o criterion, but only 5 passed 
the W2-W3 one. Two of them were identihed in the GALEX catalog, FJJ-HI and FJJ-V, 
and one, FJJ-HI, which lies close to the center of NGG 185, failed the {NUV — g)o criterion. 


The hnal dwarf galaxy, NGG 205 (d© = 824 ± 27 kpc), is projected so close to M31 
that there is some uncertainty whether a GG belongs to it or to the M31 halo. The GGs 
near the center of NGG 205, which are among the most likely candidates for membership, 
are of little use here because they he in extremely dense star helds. We therefore selected 
4 of the more remote clusters, which are also among the brightest: B009 (Hubble (H)-I), 
BOll (H-H), B020 (H-HI), and B317 (H-VHI). The B numbers are from the RBGv5, which 
provided the optical photometry. All 4 passed the criteria based on (z — lVl)o, W2-W3, and 
{NUV — g)o (asterisks in Fig. 5). 


In summary, of the 42 GGs in the dwarf Local Group galaxies discussed above, 5 (12%) 
are fainter than our tq cutoff. We could not apply our survey techniques to another 8 
clusters because of the crowding of images or the lack of g and i photometry. Twenty-nine 
GGs had adequate data and 25 (86%) passed the survey criteria. If the GGs in these dwarf 
galaxies are representative of the IGGs in the LG, then the majority of these IGGs should 
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pass our selection criteria prior to the final step of visual inspection. 


5. Results of the Local Group Survey 


As described in section 3, the color selections that are based on the SDSS, WISE, and 
GALEX measurements reduced to 183,791 objects the initial sample of millions of objects 
from the SDSS Galaxy Gatalogue. This sample is expected to contain a large fraction of the 
IGGs in the LG that are brighter than the tq limit of 19.0, if the M31 GGs are good proxies 
for them. The diagrams in Figure 2 suggest that ~ 2% of the IGGs brighter than the cutoff 
are missed by the star-galaxy separation at ~ 1100 kpc and that this decreases to ~ 1% 
at M31’s distance (783 kpc). Gonsequently, the tq cutoff appears to be the more important 
of these two limitations. It sets a lower limit on the luminosities of the IGGs that can be 
detected, which scales, of course, with their distances, and it is My < —6 at the outer 
boundary of the LG (see Fig. 2). The tests with the GGs in M31 (Figs. 3 & 4) and in the 
LG dwarf galaxies (Fig. 5) suggest that a large fraction (> 86%) that are brighter than 
the ro cutoff pass the color selection techniques. It is harder to quantify the success our 
visual inspections of the objects that pass these other criteria, particularly as a function of 
distance. The resolution of the SDSS images, typically 1'.'4 seeing with (y.'40 per pixel can 
make uncertain the distinction between compact GGs and some galaxies at the distance of 
M31 and beyond, even wi th the improvem ent of the image quality of the cutouts provided 


by SDSS in its DRIO (see 


Ahnetah 


20141 ) over previous releases. This later improvement 


would have reduced the larger misclassifications of candidate GGs that we made in our first 
survey, which was based on DR8. Our present survey retrieved 84% of the known GGs 
in M31 and M33 in the search area, which is our best estimate of the success rate of our 
methods, including the visual inspection, at ~ 70% of the radius of the LG. 


The visual examination of the SDSS cutout images of the objects that passed all 
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other criteria yielded a surprisingly small list of 17 potential new GC candidates. Five of 
them appeared enough resolved in the SDSS images that we could classify them as high 
conhdence candidates. Because they are near M31, we were able to locate them in g, r, 
or i band images taken with the MegaPrime camera of the CFHT and publicly available 
through the Canadian Astronomy Data Centre (CADC). These images are deeper and of 
higher resolution, typically 078 seeing with 0719 per pixel. With these new images, we were 
able to confirm as GCs all 5 of our high conhdence objects. Table 1 provides the positions 
of these new GCs and the data by which they were selected. The SDSS r-band thumbnail 
images of these GCs are shown in Figure 6. These faint GCs may have escaped detection 
until now because they are not within the areas covered by Hubble Space Telescope (HST) 
images and because they he at smaller Rgc than the inner boundary (~ 25 kpc) of the 
PAandAS survey, which employed MegaPrime imaging. 


To take advantage of the better resolution of the DRIO cutout images and also of the 
photometric selection criteria that are presented in this paper, we revisited the lists of high 
and lower confidence GCs provided in our first two surveys (Paper I & II). Most of the high 
conhdence objects (tables 1 of Paper I & II) passed the new selection criteria, but only three 
of the lower conhdence objects in table 2 of Paper I passed. We then visually inspected 
this reduced sample through the new cutouts of SDSS DRIO. The ones that we reclassihed 
as higher conhdence candidates could be conhrmed as GCs by our visual inspection of the 
MegaPrime imag es from the CADC archive. Some of them were later conhrmed also by the 


PAndAS search (jHnxor et ah 


20141 ). Table 2 is the updated list of our newly discovered 


and conhrmed GCs from the SDSS Galaxy Catalogue. This list is comprehensive of our 
3 surveys and counts 22 GCs, including several with Rgc > 100 kpc. It lists the positions 
of the clusters, their reddening corrected values of r and g-i (from the model magnitudes 
and extinctions listed in the SDSS), their absolute r magnitudes {Mr), their half-light radii 
{Rh), and their projected distances from M31 {Rgc)- The values of Rh were estimated from 
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the r-band light prohles provided by the SDSS, and they give rough estimates of the sizes 
of the clusters (see Paper I). 


The 12 remaining candidates that we could not confirm as GCs are listed in Table 
3, with their positions, photometric data, and angular sizes of Rh- Their SDSS r-band 
thumbnail images are shown in Figure 7. The distribution on the sky of these 12 candidate 
GCs is shown in the Hammer projection of Figure 8. There is no evidence for a concentration 
in a particular region of the sky. None of them is in close proximity to the galaxies in or 
near the LG. Two of them. Cl and C9, appear in Figure 8 to be close to the dwarf galaxies 
IC1613 and KKR-25, respectively, but in each case, the angula r separation is > 3 deg and 


many times the Rh of the galaxy given by 


McConnachie (2012). 


6. Application to the M81 Group of Galaxies 


Our survey area covers many of the galaxies of the M81 group of galaxies, including 


M81, M82, and NGC3077, which 


rave interacted with ea ch other recently and are enclosed 


m a common enve. 


Karachentsev et al 


o p of HI gas fjChvnoweth et al. 


2008 ). At distance of ~ 3.6 Mpc 


(120131 ). the M81 group is far from the LG, but nonetheless our survey 


techniques c an id entify the most luminous GCs as GC candidates. The two clusters that 
Jang et ^ ( 2012 ) identihed in the images of the HST Advanced Camera for Surveys 


(ACS), pass all our photometric criteria (hlled triangles in 
resolved into stars only on deep images with the HST (see 


Fig. 9). Since the 


VI81 GCs are 


Nantais et al. 


20111 ), our visual 


examination of the SDSS images could simply reveal if an objec t was 


or if its image resembled those of the two clusters in 


Jang et al 


obviously a galaxy 


( 2012 ). Over the area of 


the sky around the M81 group defined by 143?0 < RA < 159?0 and 64?0 < DEC < 71?0, 
we lowered the tq limit of our survey to 20.0. The SDSS footprint covers about 80% the 
region enclosed by the above coordinates and has an odd shape with the peculiar galaxy 
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M82 near its northern bonndary. Many GCs at the dist ance of the M81 gronp are 


be indistingnishable from stars in the SDSS images (see 


Perelmuter &: Racine 


199; 


ikely to 


), and 


these clusters will not be listed in the SDSS Galaxy Gatalogue. However, if the GGs in the 
M81 group resemble the M31 GGs, ~ 50% of the most luminous ones {My < —7.8) will 


be non-stellar according to the SDSS criterion (s ee Fig. 2) 


Gatalogue of the two luminous GGs identihed by 


Jang et ah 


he p resence in the Galaxy 


(1201211 is consistent with this 


expectation. The objects that we identify as GG candidates in the M81 group are plotted 
in Figure 9 (open circles) and listed in table 4, where we have also listed their photometry 
and projected distances from M81. M81-G3, M81-G4, and M81-G5 are closer to the galaxies 
M82, NGG3077, and BK6N, respectively, than to M81, and may be physically related to 
them. Two of our candidat e GGs, M81-G1 and M81-G2, are listed as 90262 and 50016 in 


Perelmuter fc Racinel (1199511 catalogue of 3774 objects within 25 arcmin of M81. Neither 


one appears to have been investigated since then. The brighter of the two, M81-G2, which 
is also closest to M81, lies within the helds of 2000s and 800s exposures in the F300W 
pass-band with the HST’s wide-held and planetary camera 2 (WFPG2). Our inspection of 


these images did not reveal any signs that M81-G2 is a galax y. The GG 1029 (INantais et ah 


20111 ). which is the most luminous one known in M81 (see 


Mavva et ah 


201311 . is very 


resolved into stars on F814W images with the HST AGS, but not obviously in a 6300s 
exposure with the WFPG2 in the F300W hlter. Gonsequently, the fact that M81-G2 is not 
clearly a GG in the F300W WFPG2 images does not rule it out as a candidate. 


7. Conclusions 

This survey for GGs in the LG has identihed 5 new GGs in the halo of M31 (see Table 
1). Their properties do not appear to be exceptional, and their discovery suggests that the 
census of GGs in M31 may be still incomplete. The searches of the SDSS Galaxy Gatalogue 
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for GCs described here and in Papers I & II have yielded a total of 22 new GCs near M31 
(see Table 2). These clusters span a wide range of Rgc (10-137 kpc), but even the most 
remote ones are likely to be members of M31’s halo. It is possible that measurements of 
their radial velocities and 3-D distances from M31 could reveal that some are not bound to 
it and are IGGs. 

Of all the objects that passed our 5 selection criteria for GGs and our visual inspection, 
only 12 are so far from M31 that they may be IGGs in the LG (see Table 3). While they 
resemble GGs in the SDSS imaging (see Fig. 7), there is still the possibility that they 
are galaxies, and we consider them to be only GG candidates. These objects need to be 
investigated in more detail to see if they are truly GGs. 

This sample of candidate IGGs in the LG is clearly incomplete because our survey has 
covered only about one-third of the sky (see Fig. 8), and we intend to expand the search in 
the southern sky. Moreover, some IGGs could have escaped our detection because they are 
hidden by obscuration near the Galactic Plane, blended with other objects, or fainter than 
our magnitude limit, which at the outer reaches of the LG excludes GGs with My > —6 
(see Fig. 2). Also, as mentioned previously, about 15% of the known GGs fail our selection 
criteria. Even with these caveats, our survey suggests that the LG does not appear to have 
a large population of IGGs, independently of their possible origin. 
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Fig. 1.— A Hammer projection of the sky with equatorial coordinates showing the area of 
our survey (gray shading), the Galactic Center (solid circle) and Plane (solid curve), M31 and 
M33 (open triangles), MW satellites (squares), M31 satellites (dots), the other LG galaxies 
(X’s), and galaxies near the LG (circles). 
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Fig. 2.— The effects of the star-galaxy separation and the magnitude limit on the surveys 
is illustrated using a sample on 168 GCs in M31. The open histogram in each diagram is 
the luminosity distribution of the whole sample. The solid histograms are the GCs that are 
non-stellar according to the criterion used by the SDSS (see text) and are brighter than the 
listed ro limit, at the assumed mean distance (<d>). Top; M31’s distance, middle: the 
outer boundary of the LG, and bottom: the distance of the M81 group. 
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Fig. 3.— The selection criteria that were used in the automated part of our survey are 
illustrated using 122 GCs in the halo of M31 (circles) and the galaxies (X’s) in a 4 deg^ 
test region of our LG survey. The dashed contours in each diagram enclose 90% of the 2-D 
density of the GGs (see text). The whole sample of GG’s is plotted in each diagram. The 
whole sample of galaxies is plotted in the top diagram, but only the ones within the rectangle 
in the top diagram are plotted in the middle diagram. Likewise, only the galaxies within the 
rectangle in the middle diagram are plotted in the bottom diagram. In the bottom diagram, 
the solid circles and the X’s represent GGs and galaxies, respectively, that have measured 
W2-W3, while the open circles and crosses are the upper limits on W2-W3 for other GGs 
and galaxies, respectively. The lines at W2-W3 = 3.1 (solid) and 4.0 (dashed) are the cuts 
for measured and upper limit values, respectively. Note that Photoz is used as a color index 
and not a measure of redshift for the GGs. 
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Fig. 4.— The sample of 122 M31 GCs (open circles) are used to illustrate the criterion that 
is based on the color formed by the GALEX NUV and the SDSS g magnitudes. The dashed 
contour encloses 90% of the density of GGs. The X’s and crosses are the galaxies that passed 
the cuts for measured and upper limit values of W2-W3 in Figure 3. The area enclosed by 
the curve and the vertical lines define the criterion used to select GG candidates. 
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Fig. 5.— The GCs in LG dwarf galaxies are plotted in the diagrams used to select GG 
candidates. The GGs in the Fornax dSph galaxy, which are labeled, are plotted with the 
WISE prohle htted magnitudes (solid circles) and the magnitude in the largest WISE aper¬ 
ture (X’s), with a dashed line connecting the two. WISE prohle htted magnitudes are used 
for GGs in NGG 6822 (solid triangles), WLM (open circle), M33 (crosses), NGG 147 (solid 
squares), NGG 185 (open triangles), and NGG 205 (asterisks). The dashed contours are the 
same as the ones plotted in Figures 3 and 4, and they enclose 90% of the density of M31 
GGs. 
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Fig. 6.— The r-band SDSS images of the newly discovered GCs in M31 (see Table 1). North 
is at the top, and East is to the left. Each image is approximately 25” by 30”. 



Fig. 7.— The r-band SDSS images of the candidate IGCs in the Local Group (see Table 3). 
North is at the top, and East is to the left. Each image is approximately 35” by 30”. 
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Fig. 8.— A Hammer projection of the sky with equatorial coordinates. The candidate IGCs 
in the LG are plotted as asterisks. The rest of the symbols are the same as Figure 1. 
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Fig. 9.— The candidate GCs in the M81 Group (open circles) and the two GGs identihed by 
Jang et al (2012) (solid triangles) are plotted in the 4 diagrams that are used to distinguish 
GG candidates from galaxies (see Figures 3 & 4). The dashed contours are they same ones 
that are plotted in Figures 3, 4, & 5, and they enclose 90% of the density of M31 GGs. 






Table 1. New M31 Globular Clusters 


Name 

R.A. 

(deg J2000) 

Decl. 

(deg J2000) 

ro 

(mag) 

{NUV - g)o 

(mag) 

{g - i)o 

(mag) 

{i-Wl)o 

(mag) 

W2 - 1V3'^ 

(mag) 

Photoz’’ 

dTZZ-06 

9.48287 

41.30932 

18.82±0.03 

2.62±0.08 

0.72±0.06 

1.11±0.15 

3.63 

0.043±0.032 

dTZZ-08 

9.63313 

39.86855 

18.22T0.01 

4.69T0.27 

1.01±0.02 

1.59±0.06 

3.28 

0.088±0.051 

dTZZ-09 

9.75310 

41.23024 

17.89±0.01 

3.89±0.07 

0.81±0.03 

1.04±0.08 

3.54 

0.119±0.046 

dTZZ-13 

11.43322 

42.65674 

18.90±0.02 

3.84±0.16 

1.04±0.03 

2.02±0.07 

3.79 

0.050±0.019 

dTZZ-14 

11.59881 

42.60960 

18.92±0.02 


1.02±0.04 

1.05±0.14 

3.62 

0.078±0.034 


Note. — (i — W1)q and W2 — W3 are on the Vega system. All other magnitudes are on the AB system. 
'^Upper limits 

'^Photoz is used as a color index, not as a measure of redshift. 



Table 2. Updated List of Newly Discovered M31 GCs in the SDSS Catalogue 


Cluster 

Name 

Previous^ 

Name 

R. A. 

(deg J2000) 

Decl. 

(deg J2000) 

ro 

(mag) 

{g - do 

(mag) 

Mr 

(mag) 

Rh 

(pc) 

Rgc 

(kpc) 

Confirmation 

dTZZ-Ol 

A 

5.14119 

36.65953 

17.44 

0.82 

-7.0 

8.5 

86 

b.c 

dTZZ-02 

B 

6.71750 

38.74947 

15.99 

0.69 

-8.5 

4.2 

54 

b,c 

dTZZ-03 

C 

7.86467 

39.53942 

17.33 

0.68 

-7.1 

4.1 

38 

b,c 

dTZZ-04 

SDSSl 

9.00774 

40.49723 

17.36 

1.08 

-7.1 

10.5 

20 

b,c 

dTZZ-05 

D 

9.03580 

39.29165 

17.24 

0.73 

-7.2 

4.1 

32 

b 

dTZZ-06 


9.48287 

41.30932 

18.82 

0.72 

-5.6 

9.7 

12 

b 

dTZZ-07 

E 

9.61483 

40.65835 

18.26 

1.01 

-6.2 

5.6 

14 

b 

dTZZ-08 


9.63313 

39.86855 

18.22 

1.01 

-6.2 

5.1 

22 

b 

dTZZ-09 


9.75310 

41.23024 

17.89 

0.81 

-6.6 

6.3 

10 

b 

dTZZ-10 

SDSS3 

9.80443 

41.70220 

18.25 

1.08 

-6.2 

5.9 

11 

b,c 

dTZZ-11 

SDSS4 

10.32496 

42.77124 

17.71 

1.14 

-6.8 

7.2 

21 

b,c 

dTZZ-12 

SDSS6 

10.61489 

39.92444 

18.51 

0.93 

-6.0 

7.1 

18 

b,c 

dTZZ-13 


11.43322 

42.65674 

18.90 

1.04 

-5.6 

4.7 

20 

b 

dTZZ-14 


11.59881 

42.60960 

18.92 

1.02 

-5.5 

6.3 

21 

b 

dTZZ-15 

SDSS8 

12.65142 

42.53047 

18.57 

1.10 

-5.9 

9.1 

26 

b,c 

dTZZ-16 

SDSS9 

13.41490 

42.58747 

17.20 

0.66 

-7.3 

7.2 

33 

b,c 

dTZZ-17 

SDSSll 

14.73495 

42.46061 

15.61 

0.70 

-8.9 

4.0 

44 

b,c 

dTZZ-18 

SDSS12 

18.19590 

42.42356 

16.88 

0.74 

-7.6 

5.0 

78 

b,c 

dTZZ-19 

SDSS15 

20.76470 

41.91971 

16.77 

0.70 

-7.7 

4.0 

103 

b,c 

dTZZ-20 

C62 

21.94838 

40.67996 

18.71 

0.76 

-5.8 

11.3 

116 

b,c 

dTZZ-21 

G 

22.20478 

47.07277 

16.98 

0.81 

-7.5 

6.1 

137 

b,d 

dTZZ-22 

SDSS16 

22.25898 

40.78570 

18.25 

0.83 

-6.2 

6.6 

119 

b,c 


^Name used in Paper I or II. 


Note. — Confirmation key: b = our visual inspection of MegaPrime images, c = Huxor et al. (2014), d = Huxor, A. (private 
communication). 







Table 3. Candidate Intergalactic Globular Clusters in the Local Group 


Name 

R.A. 

(deg J2000) 

Decl. 

(deg J2000) 

ro 

(mag) 

{NUV - g)o 

(mag) 

(g - bo 
(mag) 

(i - Wl)o 

(mag) 

(W2 - W3 

(mag) 

Photoz 

Rh 

(arcsec) 

dTZZ-COl 

13.61358 

4.18373 

18.72±0.08 

1.68±0.06 

0.31±0.12 

1.46±0.17 

2.78"" 

0.041±0.025 

2.7 

dTZZ-C02 

17.34474 

-5.91597 

18.08±0.04 


0.73±0.04 

1.67±0.06 

2.82=" 

0.072±0.075 

3.4 

dTZZ-C03 

31.37672 

6.77809 

17.84±0.01 

1.98±0.14± 

0.46±0.02 

1.31±0.06 

3.30=" 

0.034±0.010 

2.1 

dTZZ-C04 

37.649609 

46.31996 

17.25±0.01 

2.29±0.20 

0.85±0.02 

1.22±0.05 

3.02±0.33 

0.071±0.041 

3.2 

dTZZ-C05 

102.15198 

-18.38880 

17.71±0.04 


0.69±0.07 

1.19±0.10 

3.84=" 

0.100±0.057 

2.7: 

dTZZ-C06 

120.87215 

13.07730 

18.20±0.01 


0.85±0.03 

1.70±0.07 

3.95=" 

0.067±0.029 

2.2 

dTZZ-C07 

152.40421 

61.26622 

18.53±0.03 

2.01±0.07 

0.60±0.05 

0.19±0.20 

3.57=" 

0.047±0.037 

2.1 

dTZZ-C08 

238.65523 

12.92038 

18.00±0.01 

1.47±0.10 

0.37±0.03 

0.91±0.11 

3.94=" 

0.035±0.019 

2.1 

dTZZ-C09 

250.07475 

54.96822 

17.66±0.01 

3.78±0.45 

0.93±0.01 

1.41±0.03 

2.30=" 

0.074±0.036 

2.4 

dTZZ-ClO 

255.93616 

38.79727 

18.02±0.01 

2.89±0.08 

0.93±0.02 

1.90±0.04 

2.58±0.39 

0.068±0.028 

1.2 

dTZZ-Cll 

313.82498 

54.71283 

15.51±0.01 


0.85±0.02 

1.56±0.03 

1.14=" 

0.079±0.051 

1.3 

dTZZ-C12 

343.69748 

17.43929 

18.08±0.02 


0.63±0.03 

1.36±0.08 

3.90=" 

0.067±0.016 

2.8 


^Upper limit 


Note. 


Photoz is used as a color index and not a measure of redshift. 



Table 4. Candidates and Globular Clusters in the M81 Group 


Name 

R.A. 

(deg J2000) 

Decl. 

(deg J2000) 

ro 

(mag) 

(NUV - g)o 

(mag) 

{9 - i)o 

(mag) 

{i-Wl)o 

(mag) 

W2 - 

(mag) 

Photoz 

(kpc) 

M81-C1 

148.10805 

68.80711 

19.22±0.02 

1.84±0.03 

0.54±0.03 

1.89±0.10 

3.99 

0.049±0.014 

24 

M81-C2 

148.57085 

68.92275 

17.37±0.02 

3.24±0.03 

0.76±0.01 

1.46±0.04 

2.24 

0.046±0.024 

12 

M81-C3 

149.52321 

69.57946 

18.46±0.01 

3.20±0.04 

0.68±0.02 

1.57±0.08 

3.96 

0.041±0.026 

36 

M81-C4 

151.32164 

68.77572 

18.51±0.01 

3.82±0.32 

0.84±0.02 

1.08±0.11 

3.61 

0.073±0.031 

58 

M81-C5 

158.17458 

65.70965 

17.99±0.01 

2.92±0.10 

0.79±0.01 

1.74±0.04 

3.11 

0.060±0.020 

309 

GC-1'= 

148.35931 

69.52164 

18.26±0.01 

3.36±0.05 

0.76±0.02 

1.14±0.09 

3.75 

0.037±0.021 

31 

GC-2'= 

148.33411 

69.65462 

17.31±0.01 

3.29±0.03 

0.81±0.01 

1.23±0.05 

2.45 

0.067±0.030 

39 


^Upper limit 

'^Projected distance from M81. 

^Known globular cluster (Jang et al. 2012) 


Cu 


Note. 


Photoz is used as a color index and not a measure of redshift. 



